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^ ' Abstract 

■ If neutrinos are massive, 'lepton flavour eigenstates' are absent from the amplitudes of all 

Standard Model processes. Measurements of r(7r — > ev) /T{7r — > /iv) and of the MNS matrix are 
shown to exclude the existence of such states. Incoherent production of neutrinos of different 
flavour modifles the usual oscillation phase derived assuming production of a lepton flavour 
I eigenstate. The new prediction can be tested in long line-base experiments by comparison of 

'f^ disappearance' of neutrinos from pion and kaon decay. 
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^ ' In the Standard Electroweak Model, the coupling of a charged lepton, £j, of generation 

i and a massive neutrino, Uj, of generation j to the W-boson is described by the leptonic 
^ ■ charged current: 

J^(CC)''^^* = J]^,^7M(l-75)t/.,V^., (1) 

id 

where Uij is the Maki-Nakagawa-Sakata (MNS) [1] charged lepton flavour/ neutrino mass 
mixing matrix. Table 1 shows the elements of this matrix as estimated [2] from exper- 
imental measurements of atmospheric and solar neutrino oscillations. The non-diagonal 
nature of this matrix gives evidence for strong violation of generation number (or lepton 
flavour) by J^(CC)'^^*. Conservation of generation number corresponds to a diagonal 
MNS matrix with z/i = z/g, z/2 = t'^ and z/3 = Ur- This is the conventional massless neu- 
trino scenario. With massive neutrinos and a non-diagonal MNS matrix, 'z/g', 'z/^' and 
V^' do not exist as physical states. That is, they do not appear in the amplitude for 
any physical process. This is because the neutrino charged lepton couplings to the W are 
completely specifled by J^(CC)'^p* which contains only the wavefunctions of the neutrino 
mass eigenstates z/i, 1^2 and z/2. 

Indeed, the introduction of such 'lepton flavour eigenstates' as a linear superposition 
of neutrino mass eigenstates leads to predictions that are excluded by experiment. To 
show this, the pion decays: n —>■ eu and vr — > /xz/ may be considered. The invariant 
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amplitude to produce the anti-neutrino mass eigenstate Vi in association with a charged 
lepton l~ {e~ or is: 

G — 

Mti = -^Um^Vud^iil - 75)UtiipVi (2) 

where G is the Fermi constant, and rrij^ are the pion decay constant and mass re- 
spectively, and Vud is an element of the CKM [3] quark flavour/mass mixing matrix. 
Introducing a 'lepton flavour eigenstate': ip-p^ as a coherent superposition of the mass 
eigenstates ■0p.: 

i^ve = Unipvi + Ue2ipv2 + Uesipvs (3) 
gives, for the invariant amplitude of the decay process 7r~ — > i~Vf 

Me = Mf,Ua + MgUn + A^fs^^s (4) 

where Aif^ {D is for 'diagonal') is given by Eqn(2) with the replacement Uu — 1. If all 
neutrino masses are much smaller than the electron mass, it follows that: 

Mf, ^ Mg Ml ^ Ml (5) 

where Ml is the invariant amplitude calculated on the assumption of massless neutrinos. 
With two flavour {e-\i) mixing so that: 



t^el ^e.2 \ ^( COs6'i2 siu 6*12 

/ V -sin6li2 COS6I12 



(6) 



and C/e3 = t^/i3 = 0, Eqns(2)-(6) lead to the prediction: 
lip 
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'cos ^12 + sin 6112^ ^ 
cos 9i2 — sin 9i2 



(7) 



Allowing for radiative corrections [4, 5] the world average experimental value Re/f^ — 
(1.230 ± 0.004) X 10~^ [6] leads to a determination of the mixing angle Ou via the 
relation: 

l + tan^^i2\'^Qg97g^QQQ32 



1 - tan6' 



12 



The result is: 

tan2^i2 = (3.6™) x 10"^ 



or 



tan^^io < 3.0 X 10"*^at 95% CL 



These results arc evidently incompatible with the two-flavour mixing results obtained 
from the study of solar neutrino oscillations [2] where 9q ~ 9i2- 

0.22 < tan^ 9q < 0.71 (LMA solution) 

0.47 < tan^ 9q < 1.1 (LOW solution) 
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1(e) 


0.79 ±0.12 


0.57 ±0.16 


0.1 ±0.1 


2(1,) 


-0.45 ± 0.25 


0.49 ±0.28 


0.69 ±0.18 


3(r) 


0.34 ±0.29 


-0.60 ±0.23 


0.67 ±0.18 



Table 1: Values of the MNS lepton flavour/neutrino mass mixing matrix Uij as derived from solar 
and atmospheric neutrino oscillation data [2]. 



Allowing also for the coupling of the state to muons, as suggested by data on oscillations 
of atmospheric neutrinos, but still setting Ues = [7] gives, in place of Eqn(8): 



cos O12 + ^^iii 
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(cos 9i2 — sin 612) cos 623 + sin 623 



0.9976 ± 0.0032 (9) 



Assuming sin 6*12 = 1/2, and sin 6*23 = cos 6*23 = l/\/2 as suggested in Reference [7] 
(values consistent with the MNS matrix elements shown in Table 1) gives the value 2 
for the LHS of Eqn(9). It is clear, from these considerations, that the hypothesis that a 
coherent 'lepton flavour eigenstate' is produced in pion decay is experimentally excluded, 
with an enormous statistical significance, by the experimental measurements of i?e//i ^-nd 
the MNS matrix elements. 

The correct prediction of the pion decay widths, that assumes incoherent production 
of the different neutrino mass eigenstates, is: 

r(7r- ^ ru) = r(7r- ^ rz/i) + r(7r- ^ rz/2) + r(7r- ^ rus) 

~ Ml{\Un\' + \Ue2\' + \Ui3\')^Ml (10) 

where the unitarity of the MNS matrix is invoked. The predicted value of Re/^J,, given by 
setting 9i2 = in Eqn(7), is then independent of the values of the elements of the MNS 
matrix. Unlike in the case where the coherent production of a 'lepton flavour eigenstate' 
is assumed, no information about the MNS matrix can be obtained from measurements 

of Re/fi 

Although the incoherent nature of the production of the different neutrino mass eigen- 
states, as exemplified in Eqn(lO) above, was pointed out more than 20 years ago by 
Shrock [8, 9] , and the unphysical nature of coherent states of neutrinos of different mass 
was also discussed in the literature [10] the production of a coherent 'lepton flavour eigen- 
state' at a fixed time remains the basic assumption, in the literature, for the calculation 
of the phase of neutrino oscillations [7]. The phase of each mass eigenstate is then as- 
sumed to evolve with time according to the Schrodinger Equation in its own rest frame 
(the so-called 'plane- wave' approximation). The assumption that all mass eigenstates are 
produced at the same time implicitly assumes equal velocities, since there is evidently a 
unique detection event at some well defined point in space-time. Still, in the derivation 
of the phase, the neutrino velocities, as defined by the kinematical relation: v = p/E 
are assumed to be different. Thus contradictory hypotheses are made in space-time and 
in momentum space. These assumptions lead to the so-called 'standard formula' for the 
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oscillation phase [7]: 

where (i = 1, 2, 3) are the neutrino masses, L is the distance between the source pion 
and the detection event and is the neutrino momentum in pion decay to a massless 
neutrino. Units with h — c—1 are used. 

It has been shown recently [11] that the factor of two in the denominator of this equa- 
tion is a consequence of a common production time for all mass eigenstates, that follows 
from the incorrect assumption that a coherent 'lepton flavour eigenstate' is produced. As 
will be demonstrated below, allowing for the possibility of different production times, due 
to the incoherent nature of the neutrino production process, doubles the contribution of 
neutrino propagation to the oscillation phase. 

Other kinematical approximations that have been made in the literature, e.g. equal 
momenta and different energies or equal energies and different momenta, result in changes 
of only O(m^) to the phase, as compared to the exact result obtained by imposing both 
energy and momentum conservation in the neutrino production process. Thus all these 
results are equivalent at 0{rn?). 

In several previous papers written by the present author [11, 12, 13] results for the 
oscillation phase differing significantly from Eqn(ll) have been obtained. These cal- 
culations were based on the covariant Feynman Path Integral formulation of quantum 
mechanics. Here a concise re-derivation of the result for e~ production resulting from the 
interaction of neutrinos from the decay ti^ fi^u at rest will be given. For simplicity, 
the effect of the finite pion lifetime as well as that of the smearing of the physical masses 
of the TT and is neglected. As shown in Reference [12], the corresponding corrections 
to the oscillation formula are completely neghgible. The basis of the calculation is very 
simple. Each interfering path ampIitTide has a common initial state (a pion that comes to 
rest in a stopping target at laboratory time t = to) and a common final state (defined by 
the neutrino detection event). There is one such path amplitude for each neutrino mass 
eigenstate. As neither the neutrino's flavour nor its production time is observed, the cor- 
responding amplitudes must be coherently summed. Quantum mechanical coherence is a 
property of the path amplitudes (because they all have the same initial and final state) 
not of the neutrino production process. Indeed, the calculation is a direct application of 
the fundamental Eqn(7) of Feynman's original Path Integral paper [14]. 

As in the usual 'plane-wave' derivation of Eqn(ll) [7], the space-time propagators of 
the neutrino and the source pion can be simply obtained by solving the time-dependent 
Schrodinger Equation in the particle rest frame. The path amphtudes for the detection 
of an e" produced in the quasi elastic scattering process: un — > e~p by a neutrino from 
the decay at rest: tt"^ — > ii'^i/ are then: 

Ae^,{i) = AoUei eicp{-imi[TD{iy) - Ti{u)]}Uf,i exp{-im^[Ti(7r) - ro(7r)]} 

= AoC/eiC/^iexp(-iA0,) (i = l,2,3) (12) 

where: 

Ao = (e"p|T|i/n)(i/Ai^|r|7r+) (13) 

The 'reduced' transition amplitudes (e~p|T|i/n), {i'ii'^\T\7r'^) are defined by factorising 
out the MNS matrix elements, e.g.: 

(i/,A^+|T|7r+) = C/^,(i//i+|T|7r+) (14) 
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In Eqn(12), To(7r), Tj(i/, tt) and TD(i/) are the respective proper times of the tt or i/ at 
which the pion comes to rest in the stopping target, the neutrino is produced and the 
neutrino is detected. The overall propagator phase of the amplitude Agj^^i) is the Lorentz 
invariant quantity: 



A0j = mi[TD{iy) - Ti{iy)] + m^[Ti(7r) - To(7r)] 
Since the pion is at rest: 



Ti(n) — To{n) — ti — to = to — to = t^ — to — L 



2 



1 + ^ 



+ 0{nf 



where Vi is the velocity of Vi. The exact kinematical relation for the neutrinos: 

mL mL 

At = — = — = 

7 Ev p 

and Eqn(16), substituted in Eqn(15), then give: 



(15) 



(16) 



(17) 
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A0, = ^ 



2p. 



L + m^{tD -to- L) + O(m^) 



(18) 



Eqns(12) and (18) then give, for the probability of e detection distance L from the 
source: 



Pe^{L) = |^^(l)+^^(2) + ^^(3)r 



= -4|An 



U:,U;,Ue2U^2 sin^ ^ + U:,U;,Ue,U^ sm' 



+U:,U*^UeiU,, sin' ^ 



where 



and 



Am| 



L + 0(m^) 



Ami = ^'i ~ ^'j 



(19) 



(20) 



(21) 



For neutrinos produced in the decay of an arbitary particle at rest [11, 12], or the 
decay, in flight, of an arbitary ultra-relativistic particle [12], Eqn(20) generalises to: 



(22) 



where Rm = mx/ms and mx is the mass of the particle (or system of particles) recoiling 
against the neutrino produced in the decay of a source particle of mass m^. The neutrino 
is assumed to be light so that ms ^ rui. For /?-decay at rest, the oscillation phase is 
given by [11, 12]: 

Ej3 



Ami 



Pv 



2p. 



L + O(m^) 



(23) 



where Ep is the total energy release of the /3-transition. 
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The first term in the square brackets in Eqns(20) and (23) is the contribution to the 
oscillation phase of the neutrino propagator. It is factor of two larger than the similar 
contribution given by the standard formula (11). The second term in the square brackets 
in Eqns(20) and (23) is the contribution to the oscillation phase from the propagator 
of the source particle. It is opposite in sign to the neutrino contribution. It can be 
seen from Eqn(22), that, for light recoil masses such that <^ 1, the oscillation phase 
tends to vanish, strong cancellation occuring between the neutrino and source particle 
contributions. 

Examination of Eqn(19) shows that the mechanism that governs the value of Pg^t is 
interference between the path amplitudes for different neutrino flavours. A small value 
of Pg/i is not an indication of an approximate conservation of lepton flavour, but of 
strong destructive interference between the different path amplitudes, independently of 
the values of the MNS matrix elements. This is made clear by consideration of the two 
flavour mixing scenario where C/gs = and Eqn(19) simplifles, using Eqn(6), to: 

Pe^(L) = |A|2sin2 2^i2sin2^ = sin' 2^12(1 - cos 0i2) (24) 

The term — cos 0i2 in Eqn(24) originates in the interference of the path amplitudes corre- 
sponding to 1^1 and 1^2. For small values of L, e~ production is suppressed by the almost 
complete destructive interference of these amplitudes, independently of the value of 612 
i.e. of the degree of non-conservation of lepton number. The destructive nature of the 
interference is due to the minus sign multiplying sin 612 in the second row of the matrix on 
the RHS of Eqn(6). This, in turn, is a consequence of the unitarity of the MNS matrix. 

Indeed, nowhere in the description of the Vg appearence' experiment, just presented, 
do 'lepton flavour eigenstates' occur, although such an experiment is typically referred 
to [7] as 'z/^ z/g flavour oscillation'. In fact, only the mass eigenstates ui, ui and 
appear in the amplitudes of the physical processes which interfere. It is the interference 
of these amplitudes to produce the detection event that constitutes the phenomenon 
of 'neutrino oscillations'. Indeed, no temporal oscillations of 'lepton flavour' actually 
occur. Still the terms 'i/g', 'u^^ and '1/^ may still have a certain utility as mnemonics, 
even though they do not represent physical neutrino states for massive neutrinos. For 
example, it makes sense to refer to solar neutrinos, in a loose way, as a 'z/g beam ' since 
the different physical components are all created together with an electron. Similarly, 
atmospheric neutrinos are predominantly 'z/^', i.e., born together with a muon. 

The values of the correction factor C that must be applied to a value of Am'^j obtained 
using the standard formula (11), to obtain the corresponding quantity as derived from the 
path amphtude formulae of Eqns(20), (22) or (23) are presented in Table 2 for different 
neutrino sources. The correction factor is: 

0|W _ (Amf,)-*-'^ _ 1-P^ 

where Es = ms for particle decays and Es = Ep for /3-decays. 

For muon decays and purely leptonic r decays, at the mean decay neutrino energy, 
as in the flrst row of Table 2, the correction to (Am^^)**""'^ is only 16 %. The small 
recoil masses in the decays: tt — > ei/, K — > ei/ lead to very large values of C (i.e. small 
values of the oscillation phase). For the electron capture process in the last row of 
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Process 


Kinematical Condition 




11 — >• euu 
T — > ew 

T jJLVV 


/ \ 'T /OA 

{p^} = 7mi/20 


l.lD 


TT — > //i^ 




U.o72 


TT ^ ez/ 




37300 


K ^ liu 




10.4 


K ^ ev 




466718 


K TTfil' 




1.52 


T ' 


iiix = ni^ + 1)1, (C "" ' ) 


5.71 


^Na ^^-^ NA + e+ + u 


e.g. = 9.3 MeV 
for ^« Be* + e+ + u 


1.0 


e- +^ Na ^^-^ NA + iy 


e.g. e~ +^ Se — >^ Li + u 


1.0 



Table 2: Values of the correction factor C relating the oscillation phase, as calculated using the 
standard formula, to that calculated using path amplitudes for various neutrino sources. 



Table 2, which is important for solar neutrino production, the correction factor is unity. 
The same is true in /3-dccay processes when Pi^/p™"^ = 2/3, close to the mean value 
(Pi,) = {11 / 16)p'^"'^ of the spectrum of allowed /3-decays. There is then hope that mass 
differences derived from reactor neutrino experiments [15] using the formula (11) will not 
be too different from those given by the path amphtude formula (23). Actually, the phase 
should be evaluated separately for each /3-decay process contributing to the reactor flux, 
and a suitable weighted average performed. This is evidently a complex and difficult 
undertaking. In contrast the corrections for atmospheric neutrinos (predominantly from 
t: ^ fiv decay) and high energy solar neutrinos, essentially uniquely from ^B /3-decay (as 
shown in the last-but-one row of Table 2) can be applied in a straightforward manner. 

Direct experimental discrimination between the standard formula for the oscillation 
phase (11) and the path amplitude formula (22) is possible in long-baseline terrestrial 
oscillation experiments by comparing 'z/^ disappearence' as observed in neutrino beams 
produced by cither n fiu or K ^ fiu decays. Because of the relatively smaller recoil 
mass in the second process there is a strong relative suppression of the interference phase. 
For fixed neutrino momentum and source-detector distance, and considering only 'z/^ ^ z/^ 
oscillations': 

023(K) = ^023(7r) = O.O358023(7r) (26) 

Ok 

The K2K experiment with a 250 km baseline has recently found indications of 'z/^ disap- 
pearence' for neutrinos from the process tt ^ iiu oi average energy 1.3 GeV [16]. Using 
the standard formula (11) the event best fit to the data gives A.ml^ — 2.8 x 10~^ eV^ and 
023 ~ 7r/4 [16]. Inserting L — 250km in Eqn(ll) gives a phase 023 = 1.37 rad. Inserting 
these values into the two fiavour disappearence' formula [7]: 

P^^(L)~l-sin2 2^23sin2^ (27) 
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gives a suppression factor of 0.6 as compared to the no-oscillation (^23 = 0) case. Formula 
(11) predicts exactly the same phase and suppression factor for neutrinos of the same 
average momentum produced in K — /iz/ decays. However, if Eqns(22) and (26) are 
correct, (j)23{K) = 49 mrad and a reduction of only 0.06 % of the 'z/^ flux' is predicted. 
This crucial test of the basic quantum mechanics of neutrino oscillations requires the 
installation of separated low energy charged kaon beams in the JHF-SK complex [17], or 
in similar long-linebase neutrino facilities. 
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